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We present a comprehensive multiwavelength analysis of the bright, long du- 
ration gamma-ray burst GRB 070125, comprised of observations in 7-ray, X-ray, 
optical, millimeter and centimeter wavebands. Simultaneous fits to the optical 
and X-ray light curves favor a break on day 3.78, which we interpret as the jet 
break from a collimated outflow. Independent fits to optical and X-ray bands 
give similar results in the optical bands but shift the jet break to around day 
10 in the X-ray light curve. We show that for the physical parameters derived 
for GRB 070125, inverse Compton scattering effects are important throughout 
the afterglow evolution. While inverse Compton scattering does not affect radio 
and optical bands, it may be a promising candidate to delay the jet break in 
the X-ray band. Radio light curves show rapid flux variations, which are inter- 
preted as due to interstellar scintillation, and are used to derive an upper limit 
of 2.4 x 10 17 cm on the radius of the fireball in the lateral expansion phase of 
the jet. Radio light curves and spectra suggest a high synchrotron self absorp- 
tion frequency indicative of the afterglow shock wave moving in a dense medium. 
Our broadband modeling favors a constant density profile for the circumburst 
medium over a wind-like profile (R~ 2 ). However, keeping in mind the uncer- 
tainty of the parameters, it is difficult to unambiguously distinguish between the 
two density profiles. Our broadband fits suggest that GRB 070125 is a burst 
with high radiative efficiency (> 60%). 

Subject headings: gamma rays: bursts, hydrodynamics, radiation mechanisms: 
non-thermal, circumburst matter 



Introduction 



To understand the inner workings of gamma-ray burst (GRB) central engines, it is 
necessary to constrain their true energy release. If a redshift is known, the isotropic energy 
release in gamma-rays, E is0jl , is a readily measur able quantity. The "energy crisis" brought 



on by implied energy releases of > 10 erg (e.g., iKulkarni et al.lll999t iFruchter et al.lll999f ) 



was resolve d when it was realized that GRB blast waves are collimated with ope ning angles 
6j of 2-30° jRhoadslll999l : ISari et al.lll999l : lstanek et al.lll999l : lHarrison et aDl999h . Thus the 
beaming-corrected gamma-ray energy release, E 7 , is smaller than -Ei SO , 7 by a factor ~ $f/2. 

Whatever energy is not released in the prompt emission powers a relativistic blast wave 
that plows into the circumburst medium (CBM). Measuring the kinetic energy of this outflow, 
Ek, is challenging because only a fraction of this en ergy is radiated. Several methods have 
been used in this endeavor (see iBerger et al.ll2003bl ). The simplest method, using only the 
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X-r ay afterglow l i ght cur ve, was suggested by iKumarl (120001 ). iFreedman fc Waxmanl (120011 ). 
and iBerger et al.l (j2003al ). Provided that the X-rays are predominately synchrotron and the 
electrons radiate in the cooling regime, this method yields the energy per unit solid angle 
and the fraction of the shock energy carried by electrons e e . For those afterglows with high 
quality multiwavelength data sets, a fit of the data can be made using models that describe 
the dynamics of jet/circumburst int eraction and that calculate the expected sy nchrotron 
and inverse Compton emission (e.g., iPanaitescu fc Kumar! l200ll ; lYost et al.ll2003l ). Finally, 
for bright, long-lived afterglows, can be measured more robustly using the late-time 
radio light curve. Months or in some cases years after the burst, the blast wave becomes 
sub-relativistic ( i.e. Sedoy self-similar evolution) and the outflow i s expected to be quasi- 



spherical (e.g., iFrail et al.ll2000l : IBerger et al.l 12004 ; iFrail et all 120051 ). 



Early studies, based on these methods, suggested that the energy release of l ong-duration 



GRBs lies within a narrow range with a mean of ~ 10 erg (e.g. IFrail et al.ll200ll ; IBerger et al. 



2003al ). There are several problems with this simple picture, however, and the true situation 
is likely more complicated. For example, energy input into the blast wave now appears to 
extend beyond the prompt emission and into the afterglow phase. This could be due either 
to a l ong-lived central en gine or slow-moving ejecta that take a long time to reach the shock 
front Jzhang et al.ll2006h . 



The simple pictures of collimated explosions has been challenged by Swift detected 
GRBs. Kocevski & Butler (2007) give an excellent summary of the difficulties and their 
implications for collimation-corrected energy release E~. If (for whatever reason) the Swift 
sample of bursts has no jet breaks then we are faced with a growing number of GRBs that 
(by virtue of inferred is otropic energy release) are hyper-energetic i.e. well in excess of the 
canonical 10 51 erg (e.g.. iBurrows fc Racusinll2007bl ; ICenko et al.ll2006bl ; IFrail et al.ll2006l ). If 
accepted, these events would provide stringent tests of existing progenitor models. 

In this paper, we present multiwavelength observations of GRB 070125. Its bright af- 
terglow has allowed us to follow the GRB until day 350 and obtain the most extensive radio 
data in the Swift era, coupled with well-sampled X-ray and optical light curves indicative of 
a jet break. Toget her with the well characterized prompt emission extending beyond 1 MeV 
dBellm et al]l2007h . GRB 070125 is truly a rare event. 



In f|2]we provide details of observations for GRB 070125. In $3] we describe our results, 
and find evidence for a jet break in the optical and X-ray data. The radio spectra show 
evidence for evolution from an optically thick to an optically thin phase, while the radio 
light curves show short timescale variability that we ascribe to scintillation. We carry out 
a straightforward analytic modeling and derive some physical parameters of the shock and 
the circumburst medium. Finally, we carry out a detailed model fit of the entire broadband 
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dataset and summarize our results in SQl We discuss the energetics, environment, relevance 
of reverse shock emission, efficiency of the GRB and the inverse Compton scattering effects 
for GRB 070125 in §5J The main conclusions are listed in £0 



2. Observations 

In this section, we present multiwavelength observations of GRB 070125. We supplement 
our data set with the measurements reported in literature (mostly notices from the Gamma- 
ray Burst Circulars Network^). 



2.1. Gamma-Ray Observations 



GRB 070125 was discover ed by the Inter Pla netary Network (IPN) of GRB detectors at 
07:20:45 UT 25 January 2007 Iflurlev et alJl2007t ). Mars Odyssey (HEND and GRS), Suzaku 
(WAM), INTEGRAL (SPI-ACS), RHESSI and Konus-Wind, all observed this intense, ~ 70 
s long, event (Fig. [T]). The burst was not in the field of view of Swift. Four minute s later 
the burst fell into the Burst Alert Telescope (BAT) field of view (IHurley et al.l 120071 ). The 
Swift BAT position was consistent with the position of the GRB as triangulated by IPN. 
The small 4' error circle from the BAT enabled follow-up observations to identify the bright 
X-ray (§2J2]) and optical (|g3J afterglow. 

Because of the high-energy (up to 1 MeV) coverage provided by RHESSI, Konus-Wind, 
and Suzaku, tight constraints can be made on the prompt emission spectrum of GRB 070125. 
In particular, the gamma-ray fluence was much more accurately measured than a typical Swift 
GRB (with coverage of only ~ 300 keV). A detai l ed an alysis of the high-energy properties of 
GRB 070125 has bee n performed by iBellm et al.l (120071 ). who found that the spectrum is well 
fit by a Band model (IBand et al.lll993l ) with a peak energy E p = 430 keV and the power-law 
photon indices (as defined in Band model) of —1.14 and —2.11 below and above the peak 



energy, respectively. The resulting fluence in the 20 keV - 10 MeV band was 1.7 x 10 

-2 



erg 



cm 



1 http:/ /gen. gsfc.nasa.gov/gcn3^archive. html 
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2.2. X-ray observations 



The X-ray Telescope (XRT; iBurrows et al.ll2005af) on board Swift began observing the 
field of GRB 070125 46 ks after the burst trigger (IRacusin fc Veterel l2007h . Due to the 
relatively low source flux, observations were conducted exclusively in Photon Counting (PC) 
mode, and so the effects of photon pile-up were negligible. During the first five orbits, the 
X-ray spectrum is well-fit by a power-law model with a photon ind ex T = 2.1 ± 0.3 (AT (7) oc 
7 _r G?7, where iV is the photon flux density, iRacusin fc Veterell2007l ). The i nferred absorption 



is con sistent with the Galactic absorption value of Nr ~ 8 x 10 cm ( IRacusin fc Vetere 



20071 ). The XRT continued to monitor the X-ray afterglow of GRB 070125 over the course 
of the next two weeks, until the source faded below the XRT sensitivity threshold. 



Motivated by the proposed lack of X-ray jet breaks observed in the Swift era (IBurrows &: Racusin 



2007bl ). we obtained Director's Discretionary time on the Chandra X-ray Observatory to ob- 
serve the X-ray afterglow of GRB 070125 at very late times. We obtained a 30 k s exposure 
using the Advanced CCD Imaging Spectrometer on board (IGarmire et al.l 120031 ) beginning 
at 21:28 UT 2007 March 5 (~ 40 d after the burst trigger). No source was detected at the 
position of the X-ray afterglow at this time. Formally, using a circular aperture with a one- 
arcsecond diameter, we detect 0.9 ± 5.0 photons in the energy range from 0.3 — lOkeV at 
the location of GRB 070125 (a 5 — o upper limit). 

The log of X-ray observations and the measured fluxes are s ummarized in Tab le [TJ 
The Swift XRT light curve is obtained from the on-line repositorj@ ( Evans et al. 2007 ). We 
converted counts to flux using spectral parameters derived from the first five XRT orbits, 
ranging from 47 — 389 ks. There is no evidence for spectral evolutio n in the X-ray light curve, 
implying the derived conversion factor is applicable at all times ( IRacusin et al.l 120071 ). We 
then converted the 0.3 — 10.0 keV flux to a flux density (F u oc v~@, where (3 = V — 1 = 1.1) 
at E = 1.486 keV (vq = 3.594 x 10 17 Hz). This value was so chosen that the flux was equally 
divided below and above this cutoff. 



2.3. Optical observations 



In response to the IPN-Swift localization, we bega n observing the fiel d of GRB 070125 
with the automated Palomar 60-inch telescope (P60; ICenko et al.l l2006al ) on the night of 
2007 January 26. Inspection of the first images revealed a bright, stationary source at 
a = 07 h 51 m 17"75, S = +31°09'04"2 (J2000.0), not present in the Sloan Digital Sky Survey 



2 http: / /www. Swift. ac.uk/xrt_curves 
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(SDSS) images of the field (lAdelman-McCarthv et al.l 120070. This position was promptly 
reported as the afterglow of GRB 070125 (jCenko &: Fox! 120071 ). allowing several groups to 
obtain spectroscopy of the afterglow while still quite bright (R ~ 18 mag). 

The optical afterglow of GRB 070125 was observed by a variety of facilities worldwide, 
as it turned out to be one of the brightest optical afterglows ever detected for a GRB (Updike 
et al, in prep.) We continued to monitor the afterglow of GRB 070125 with the P60 in the 
Kron R and Sloan i' filters for the following four nights, until the afterglow faded below 
our sensitivity lim it. All P60 data were reduced using our custom software pipeline (see 
Cenko et al.ll2006al for details) using iraf§ routines. 



In addition to our P60 monitoring, we obtained three epochs of late -time optical pho- 
tometry with the Gemini Multi-Object Spectrograph and Imager (GMOS; iHook et al.ll2003l ) 
mounted on the 8-m Gemini North telescope. All three epochs consisted of either single r' 
or i' exposures obtained as acquisition images for spectroscopic observations. All GMOS 
images were reduced using the IRAF gemini package. 

Finally, we obtained a single epoch of simultaneous g- and i?-band imagi ng on the 



night of 2007 February 16 with the Low-Resolution Imaging Spectrometer (LRIS: lOke et al. 



19951 ) mounted on the 10-m Keck Telescope. Individual images were bias-subtracted and 
flat-fielded using standard IRAF routines. Co-addition was performed using Swarpfl Pho- 
tometric calibration for all of our opti cal imaging was p erformed relative to the SDSS, with 
empirical filter transformations from iJordi et al.l (120061 ) applied where necessary. Typical 
RMS photometric uncertainties were ~ 0.05 mag in all of our images. 



Garnavich et al.l (120071 ) imaged the position of the GRB 070125 afterglow with the LBC- 



blue CCD camera and 8.4-m SX mirror on the LBT on 2007 February 21.1 (UT) in r-band. 
A faint source is detected at the position of the afterglow with brightness r=26.3 ± 0.3 mag. 
Since the source may be contaminated by the host galaxy, this observation represents an 
upper-limit on the magnitude of the afterglow 26.8 days after the GRB. 



T o convert magnitudes to flux densities, we used zeropoint measurements from lFukugita et al 
dl995]). We have incorpo r ated a modest amou nt of Galactic extinction (E(B — V) = 0.052; 
Dickey fc Lockmanlll990l ; ISchlegel et al.lll998l ) into these results. The results of our optical 
monitoring of GRB 070125, together with measurements reported by other observatories via 
the GCN, are shown in Table [2j 



3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associ- 
ation for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

4 http: / /terapix. iap.fr 
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2.4. Millimeter & Sub-Millimeter observations 



We obtained director's discretionary time for observations in the 1.2-mm band (250 
GHz) at the Max-Planck Millimeter Bolometer Array (MAMBO), installed at the IRAM 30 
m telescope on Pico Veleta, Spain. We used the MAMBO-2 version with 117 channels. The 
bandwidth used was 210 — 290 GHz half power. Our first observations took place on 2007 
January 30 and we detected the afterglow of GRB 070125 at a flux density of 3.14 ± 0.59 
mJy We monitored the afterglow regularly until it had dropped below the instrumental 
sensitivity. 

Observations were also obtained at 95 GHz using the Combined Array for Research in 
Millimeter-wave Astronomy (CARMA jl in good weather on 2007 February 5, 9, 12, and 
18. Individual observations were between 6.7 — 8 hours in length. System temperatures 
ranged from 200 to 400 K, scaled to outside the atmosphere. A single linear polarization was 
received in a bandwidth of 1.2 — 1.5 GHz per sideband, depending on the antenna. CARMA 
was in the C configuration, with baselines 26 — 370 m. 

The bright quasar 0748+240 was observed as a phase calibrator at 20-minute intervals. 
A fainter source nearer GRB 070125, 0741+312 (2° distant), was observed for 30 seconds 
after each observation of the GRB. The flux density of this source when self-calibrated was 
compared to the flux density when calibrated with 0748+240. The ratio of these measure- 
ments provides an estimate of the atmospheric decorrelation affecting the observations. The 
measured flux densities of GRB 070125 were corrected by this factor, which ranged from 9 to 
14%. The absolute flux density scale was derived from observations of Uranus on February 2 
and 19, and transferred by reference to the quasar 3C84, which was observed on those dates 
and near the beginning of each observation of GRB 070125. 



Images were made for each observation using the MIRIAD software package (ISault. Teuben. fc Wright 



19951 ). The quoted uncertainty in the measurement is dominated by the sensitivity of the 
map, but also includes an estimate for the flux density scale uncertainty of 15%. We de- 
tected the afterglow in the first three observations. The flux density of the mean observation 
time 2007 February 5, 0700 UT was 2.3 mJy The log and flux densities of the MAMBO and 
CARMA observations can be found in Table [3j 



5 Support for CARMA construction was derived from the Gordon and Betty Moore Foundation, the Ken- 
neth T. and Eileen L. Norris Foundation, the Associates of the California Institute of Technology, the states 
of California, Illinois, and Maryland, and the National Science Foundation. Ongoing CARMA development 
and operations are supported by the National Science Foundation under a cooperative agreement, and by 
the CARMA partner universities. 
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2.5. Centimeter band observations 

The earliest measurement of the radio flux dens ity was taken from the Westerbork 



Synthesis Radio Telescope (WSRT) in 5 GHz band by Ivan der Horstl (120071 )). just 1.5 days 
after the burst. A 2— a upper limit on the 5 GHz flux density of F u < 174 /iJy was obtained. 
Shortly thereafter, we triggered Very Large Array (VLA) observations of the field. Our 
first measurement at 8.46 GHz, four days aft er the explosion, resulted in a strong detection 



{F v = 360 ± 42/zJy; IChandra fc Fraill 120071 ). Encouraged by this detection, we triggered 
observations with the Giant Meterwave Radio Telescope (GMRT) in the 610 MHz band 
on 2007 January 31. However, we di d not detect the afte rglow with the GMRT to a 2— a 



limiting flux density of F v < 300yuJy (IChandra et al.l 120071 ). 



We continued the followup observations of GRB 070125 with the VLA from 2007 Jan- 
uary 29th until 2007 February 2nd at 8.46 GHz and 4.86 GHz. After 2007 February 5, we 
undertook observations in the 1.46, 4.86, 8.46, 14.94 and 22.5 GHz bands. We followed the 
GRB until day 342 since explosion. 

Each observation at a single frequency was from 30 minutes to one hour duration. The 
bright radio quasar 3C48 (0137+331) was used as a flux calibrator. We used phase calibra- 
tors 0745+317 and 0741+312 to track the instrumental and/or atmospheric gain and phase 
variations, as well as to monitor the quality and the sensitivity of the data. A bandwidth 
of 2 x 50 MHz was used for all the observations. The data were analyzed using standard 
data reduction procedures within the Astronomical Image Processing System (AIPS). During 
initial observations, the VLA was in "C^D" and "D" array configuration for early observa- 
tions; because of the short baselines, there were many sources in the field of view at lower 
frequencies. For such cases at 4.86 and 1.43 GHz, we did 3-D cleaning with several rounds 
of self-calibration inside AIPS. 

On 2007 February 7, 8 and 14, we made longer observations for durations of ~ 5.5, ~ 4, 
and ~ 8 hours, respectively, in order to search for variability due to Interstellar Scintillation 
(ISS) (see §3.2p . We combined every 20 minutes of the data and imaged the afterglow field 
of view to measure short timescale scintillations. The results of our radio monitoring of the 
afterglow of GRB 070125 can be found in Table [3j 



3. Results 



In this section, we carry out a simple analysis with minimum model assumptions to 
derive robust conclusions. We find strong evidence for a jet break in the combined fitting of 
the X-ray and the optical light curves ( §3.ip . In §3.21 we use the radio scintillation data to 
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constrain the size of the emitting region. Together, the jet break time and the inferred radius 
allow us to constrain the prompt energy release ( §3.31) and the density of the circumburst 
medium ( §3.41) . In the next section (fgj), we undertake a detailed analysis utilizing the full 
machinery of afterglow models. 



3.1. Break in the Light Curve 

We performed a joint fit on the R, i', and X-ray light curves of GRB 070125 using both 
a single power-law (spl) and a broken power-law (bpl) model. The results of these two fits 
are shown in Figure [2] and Table HI We also perform the spl and the bpl fits letting the 
X-ray and optical pre-break indices vary independently. The results are consistent with the 
previous fits in which we constrained both the indices to be the same. Altogether, the broken 
power-law model, indicating a jet break at t- s = 3.8 d, is strongly favored (x 2 (bpl) = 1-30 for 
97 d.o.f. vs. x r 2 (spl) = 2.23 for 99 d.o.f.). 

We note, however, that the distinction between the two models is significantly more 
pronounced in the optical (x 2 (-R— band, bpl) = 1.54 for 52 d.o.f. vs. band, spl) = 

2.22 for 52 d.o.f.; and x r V-band, bpl) = 0.77 for 25 d.o.f. vs. y 2 A i ' -band , spl) = 3.43 



for 52 d.o.f.) than the X-rays. As first noted by iBurrows fe Racusinl (l2007al ). without the 
late-time Chandra data, the X-ray light curve cannot distinguish between the spl and bpl 
models. Formally, we find the single power-law model is actually favored in the X-rays 
(x r 2 (X-ray,spl) = 0.88 for 22 d.o.f. vs. x?(X - ray, bpl) = 1.32 for 20 d.o.f.). This could 
perhaps be due to the denser sampling at early times. 

We also perform independent optical and X-ray fitting. The optical fits are consistent 
with our joint fits. However, when we fit the X-ray light curve independently of the optical, 
the jet break appears much later, t- } ~ 9d (Fig. [2]). While not formally required by our 
simple analysis, we consider the possibility that this break may in fact be chromatic. We 
discuss this further in SJ5J 



3.2. Scintillation and fireball size 

As can be seen in Fig. [3] and Table [31 there are significant day-to-day deviations in the 
low frequency radio light curves. These modulations are likely due to scintillation caused by 
interstellar propagation effects. 

We obtained long-duration observations of the afterglow of GRB 070125 at 8.5 GHz on 
three separate occasions: February 7 (5.5 hour duration), February 8 (4 hour duration), 
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and February 14 (8 hour duration). The data were split into 20 minute blocks and imaged 
in order to extract information on the fast variability of the source. GRB 070125 exhibits 
flux density variations with a significance exceeding 99.8% on each of the three epochs (see 
Fig. H]). The reduced \ 2 values for the hypothesis of no variability are 2.31 (with 17 degrees 
of freedom), 4.42 (10 d.o.f.) and 2.84 (23 d.o.f.) for the three dates, respectively. 

We used intensity structure functions to determine the variability timescale on each 
date. While the short baseline hindered a definitive determination, we tentatively identify 
breaks with 20 — 30% accuracy at AT « 6 x 10 3 , 7 x 10 3 and 9 x 10 3 s, in the three structure 
functions. The breaks are marginally significant in the data from 7 February and 8 February, 
but not so prominent in the 14 February images. This can be explained by quenching of the 
scintillation as the source expands at late times. These measurements are congruent with the 
timescales of the peaks and troughs apparent in the corresponding centimeter-wavelength 
lightcurves (Fig. [3]). 

The interpretation of the variability depends on whether the scintillation occurs in the 
weak or strong regime. Strong scintillations require that the so-called Fried parameter (co- 
herence length scale, s^) be smaller than the Fresnel size (rp). Strong chromatic scintillations 
i s possible only for sources smaller than X/sd- These lead to the following two conditions 



f lGoodman!ll997h 



SM \ 6/17 /A ^ 5/1? 



- <13 - 4 l lO—m- 2 o/3 kpc j {^) GHz = z/ ss (1) 

( v V 11/5 f SM V /5 / c 

where SM is the scattering measure and D SCT is the effec tive distance to the s cattering 



material, which is essentially the Galactic ionized medium (ICordes &: Lazioll2002l ). To de- 



termine this, we first estimate th e scatt ering distance and the scattering measure using 



the formulation of ICordes Laziol (120021 ). Given the Galactic coordinates of GRB 070125, 
(l,b) = (189.4,25.6), the expected SM is 10" 319 m~ 20 / 3 kpc and the effective distance to the 
scattering material is D scr = 0.84 kpc. For these parameters, the critical frequency obtained 
from Eq. [TJ is 12.7 GHz. These parameters also satisfy the condition of Eq. Hence, the 
GRB is likely to be in the strong scintillation regime. However, given the paucity of our 
knowledge of the distribution of scattering material off the Galactic plane, these estimates 
should be taken with caution. 



Strong scintillation can be diffractive as well refractive in nature. Diffractive scintillation 
of a point-like source is characterized by variations with a modulation index (m p ) close to 
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unity, where the modulation index is a measure of how much the modulated flux varies around 
its intri nsic flux density level. The size limit for a source to exhibit diffractive scintillation 
(6 X ) is JGoodmanlll997r i: 



9 1 = 1.2 



8.46 GHz 



6 / 5 (D R 



kpc 



SM 



10 -3.19 m -20/3 kpc 



-3/5 



/ias. 



(3) 



For actual size of the source, 9 S , 9 S < 9x- Diffractive scintillation occurs when 9 a > 9x, which 
is indeed the case for GRB 070125 using the parameters discussed above. 

We now atte mpt to limit the size of the emitting regio n using the formulation of diffrac- 
tive scintillation (jWalkerlll998l . l200ll ; ICordes fc Lazioll2002l ). but with more robust estimates 
of D SCT and SM. The timescale of diffractive scintillation can be expressed as: 



tdiff = 5950 



v \ 6 / 5 / SM 
8.46 GHz J Vl0- 319 m- 2 °/ 3 kpc 



-3/5 



^ISS 



30 kms 



-i 



(4) 



where t>iss is the speed of the scattering material transverse t o the line of sigh t. The decor- 
relation bandwidth (Az/d c ) for diffractive scintillation is then (]Goodmanlll997l ): 



Az/ dc = 1.55 



v 



\ 2 2/5 ( D 

8.46 GHz J 



kpc 



SM 



10 -3.19 m -20/3 kpc 



-6/5 



GHz. 



(5) 



Using the variability timescale of AT ~ 7 x 10 3 s measured on 8 February, along with 
t>iss = 30 kms -1 , we infer a scattering measure of SM = 0.76 x 10~ 319 m~ 20 ' 3 kpc. In Eq. [5] 
we take the decorrelation bandwidth to be roughly half of the frequency of observation, 
i .e. A^Hr ~ (1/2 ) x z/, because at the boundary of strong and weak scintillation Au^ « v 
( IGoodmanlll997l ). This yields a distance to the scattering screen of D scr = 0.51 kpc. 



Knowing these two parameters, we can put an upper limit on the angular size of the 
emitting region of # src = 2.8 ±0.5 //as. Here t he error in angular size corresponds to 20% error 
in the determination of AT. At z = 1.547 ( Cenko et al. 2008 ). the angular size translated 
to a linear size of (4.7 ± 0.8) x 10 17 cm or a radius of 2.4 x 10 17 cm. 

Variability due to refractive scintillation is expected at late stages, even after the source 
has expanded sufficiently Refractive scintillations are expected to be broadband in nature 
and start dominating once diffractive scintillations are quenched. The source exhibited 
refractive scintillation during the two months subsequent to 2007 February 8, as can be seen 
from the first 70 days of the 8.46 GHz data plotted in Figure [31 



3 All the calculations were done with Hq = 71, fl m = 0.27 and Vl vac = 0.73. 
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In a standard GRB afterglow model, the jet starts to spread sideways after the jet break 

( jPiranlll999l . 120051 ; lMeszarosll2006l ) . During this stage, the fireball size remain constant. Once 
the jet has become spherical, it reaches the non-relativistic regime. In this regime, the equa- 
tions of motion follow self-similar Sedov- Taylor solutions. Thus, this estimate represents the 
size of the fireball in the post jet break regime until the expansion becomes sub-relativistic, 
which occurs between days 30 and 50. 



3.3. Energetics 



The isotropic gamma ray energy of a GRB can be written as: 

2 
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where &l is the luminosity distance of the GRB and / 7 is total fluence. The fluence in the 
20 keV - 10 MeV energy range is 1.74 x 10~ 4 erg cm -2 (IGolenetskii et al.l 120071 ; iBellm et al. 



20071 ). yielding E ltiso = 1.06 x 10 54 erg. This is one of the largest isotropic energy releases 
(top 1%) ever reported for a GRB (jAmatill2006l ). 



To determine the true prompt energy release, this isotropic value needs to be corrected 
for collimation. The combined fit to the optical and X-ray data gave a jet break at £j w 3.78 
day ( §3.ip . The collimation correction depends on the density profile of the circumburst 
medium. We derive corrections for a uniform density (ISM; n = constant) medium and a 
wind-like (wind; n = Ar~ 2 ; A = 3 x 10 35 At cm -1 ) medium (jWu et al.ll2005l. and references 



therein). The collimation angle for a radiative afterglow can be written as (jSari et al.l 11998 
Frail et al.lfioOll : Il4 fc Chevalier] l2003h : 
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where tj is the break in the optical light curve in days, and E52 is the isotropic kinetic energy 
of the fireball in units of 10 52 erg. The term T is the initial Lorentz factor of the fireball. 

Similarly, the size of the spherical fir eball can be expressed in the ISM and wind media 
as JSari et alill9~98l : lli fc Chevalieri 12003 ^ : 



#(ISM) = 1.3 x 10 



17 
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«(W M ) = 0,7 x 10" (f ) 2 ' 3 (I^y-* (i±i)" /3 C on, (9) 
which after jet break in GRB 070125 translates to: 

fl ,ISM,^,x 1 0"(^) I/3 (^)" ,S (lL) ,/3 c m (10) 

TO H0,7xl0"(f)(^)"(^) cc (11) 

Using the emission radius derived from scintillation studies (R < 2.4 x 10 17 cm), we find an 
opening angle of 9 < 0.25 rad (14°) for the ISM model and 9 < 0.23 rad (13°) for the wind 
model. The Lorentz factor of the shocked ejecta at the time of jet break is 7(tj e t) — 1/6* = 4 
in the ISM model and 5 in the wind model. 



3.4. Circumburst density 

From Eq. (jTOl) and (TTTT) . the circumstellar density can be written in terms of the kinetic 
energy of the afterglow as n > -E52/O.II cm -3 and A* > E 52 /1.5A for the ISM and the wind 
models, respectively. Let rj 1 be the efficiency factor for converting the fireball ener gy into the 



radia tion energy, i.e. ^ 7 = E^/(E^ + E K ). For an empirical value of ?7 7 = 0.35 (IFrail et al. 
20031 ). the number density of GRB 070125 is n « 50 cm -3 in the ISM model and A* « 2.5 in 



the wind model. This value is quite high, even for GRBs, and indicates that the afterglow 
of GRB 070125 is expanding into a dense medium. 

A natural consequence of a high circumburst density is a high synchrotron self-absorption 
frequency. To this end, we plot broadband radio spectra (Figure [Sj). To improve the signal- 
to- noise ratio, we binned the spectra into 4 groups (t = 6 — 17 d, 18 — 35 d, 56 — 86 d, and 
160 — 200 d). The division was initially done on the basis of similarly-looking spectra, but 
later in this section we justify it by demonstrating that the spectral evolution has a weak 
time dependence. As can be seen, there is a clear turnover in the spectra in the first three 
epochs. Moreover, there is some indication that the turnover frequency evolves to lower fre- 
quencies with time. At the last epoch, no turnover frequency is discernible and the spectrum 
is inverted. 

This spectral behavior has been seen in many previously well-studied GRB afterglows. 
We interpret this behavior in terms of the evolution of the afterglow from an optically thick 
to optically thin phase, parametrized by a single unknown, the synchrotron self-absorption 
frequency v a . We can measure the value of v a from the radio spectra using a very simple 
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formulation given below: 
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The above relation is a broken power-law with a break at u a . We used the following smooth 
approximation of eq. (|T2|) : 



F = F 

v A max 



1 + 



2-^ 



(13) 



We fit this function to the first three radio spectra and obtain the following values of syn- 
chrotron self-absorption frequency: 



12.25lJ;^ GHz in the range 6-17 days, 
v a = 11.22±J;g? GHz in the range 18 - 35 days 
v a = 7.49±g;|| GHz in the range 56 - 86 days. 



Our approximation of a constant v a within each epoch is justified because of the slow 
evolution of v a . The best fit time dependence to the values given above is u a oc £~°- 24±0 - 05 , 
which agre e s wel l with the ti me d ependence predicted in the fireball model {u a oc t~ a2 , 



Frail et al.l (120031 ): iMeszarod (120061 )). The spectrum in the final epoch was well into the 



optically thin phase so we could not determine v a at this epoch. 

We also plot the multiwaveband spectra on day 10.7 and day 23.4 (Fig. [6]), the two 
epochs at which we had observations in all the bands simultaneously. This shows radio data 
points to be on the optically thick part of the spectra. The afterglow peaks at 3-mm band 
and the optical and X-ray data fall in the optically thin regime. This gives a rough estimate 
of the break frequencies, v my corresponding to the minimum electron Lorentz factor, and u c , 
the cooling frequency. We also plot vF v , a measure of energy, against v. The curve peaks 
at the cooling frequency v c (Fig. [6]). 



4. Broadband Modeling 

In the previous section, we used simple analytical techniques to estimate four funda- 
mental physical properties of the explosion: the opening angle (6 fa 0.23 - 0.25 radian), the 
size of the emitting region (R fa (2.4 x 10 17 cm), the collimation-corrected prompt energy 
release (J5 7 ~ 3 x 10 52 erg), and the circumburst density {n ~ 50cm -3 ; A* fa 2.5). Here 
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we combine all our observations of the GRB 070125 afterglow in an attempt to derive a 
comprehensive model of the entire afterglow evolution. Our modeling software assumes a 
standard synchrotron forward shock formulation, with possible contributions from inverse 
Compton (IC) emission and radiative losses also included. This model also includes scintilla- 
tion uncertainti es, and henc e giv es realistic e stima tes of various parameters. Further details 
can be found in lYostl J2004h and lYost et aD J2003h . 

As in §3.11 we ignore all data before t = 1 d due to the possibility of late-time energy 
injection. Based on our results in §3.2[ we have incorporated scintillat ion effect s into the radio 
regime. We used an LMC-like extinction model for the optical data (|Peilll992l ); however, the 
extremely small host contribution to the extinction makes the effects of differing extinction 
laws negligible. 



4.1. Wind Model 

Results of the best fit wind model parameters are tabulated in Table O In terms of 
X 2 and the model-fit statistic]^, the wind model does a slightly better job than the ISM 
model. However, the resulting best-fit parameters for the wind model are either unphysical 
(i.e. e e , the fractional energy imparted to electrons in the shock, approaching unity) or quite 
different from the values we derived in §31 The extremely small isotropic afterglow kinetic 
energy (E52 ~ 0.3) compared to the 7-ray isotropic energy is also troubling (the same is true 
in the ISM case, but to a lesser extent). We also notice that the wind model is less stable to 
small changes in the parameter space. 

We fit the wind model with and without the inverse Compton effects. The model without 
inverse Compton effect gives even more unphysical values with many of the parameters 
asymptotically reaching very high values. The magnetic field fraction required reaches 100% 
in this model. 

Apparently allowi ng e P to be a free p aram e ter is problematic. Microphysics evolution 



has been considered by lYost et al.l (120031 ); lYostl (120041 ) . It makes everything unconstrained. 



We fixed e e to be 0.4 and obtained a good fit (see Table [5]). This exercise demonstrates that 
microphysical parameters are not constrained by our observations, at least for wind model. 



7 — ln(-P) = 0.5 x (x 2 + 2£ Zn(cri))+constant, here P is the probability function and ui is standard deviation 
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4.2. ISM Model 

The best-fit results for the ISM model are also provided in Table |5j The ISM model 
gives values of various parameters closer to the ones obtained from our simple analytical 
models. The jet break time, density scale, and collimation angle are in good agreement with 
our previous results. 

The energy quoted in Table[5]is the isotropic blastwave energy at the time when v c = u m , 
i.e. at the time of the transition from fast cooling to slow cooling (t ~ 8d in our model). 
This isotropic kinetic energy is much smaller than the isotropic gamma-ray energy obtained 
from the Konus-Wind/RHESSI fluence ( §3.31) . This may indicate that either there are high 
radiative losses at early times or the prompt emission is rather efficient with an extremely 
high value of rj 7 . 



4.3. Broadband Model Results and Interpretation 

We plot the results of our broadband modeling in Figures [7J El and M It is difficult to 
differentiate between the wind and the ISM models purely on the basis of these plots. Both 
models represent the optical data fairly well at early times (Fig. [7]). At late times (t > 4d), 
both models over-predict the .R-band flux, though with significantly less discrepancy in the 
constant density medium. IC effects are negligible in optical bands. 

In the radio bands, none of the models fit particularly well, especially at early times 
(< 15 days). This is likely caused partly by the diffractive scintillation ( §3.2)) . IC scattering 
has no influence in this band. The most puzzling behavior is revea led in the 4.8 GHz ba nd. 
An early detection at t ~ 4d, both in our data and the WSRT Ivan der Horstl (120071 ). is 
followed by almost two weeks of non-detections. Summing all these non-detections, we can 
put very strict limits on the 4.8 GHz flux at this time: /„ < 71 /iJy. The reason for this 
drop in flux is unclear, for it cannot be due to scintillation. We consider the possibility that 
the early detection is due to the reverse shock emission, in §5j 

In Figure [SJ we plot the broadband spectra from radio to X-ray at various times of the 
afterglow evolution. They are represented well with our models. Both the wind and ISM 
environments seem to do an equally good job. 

We plot the X-ray light curve of GRB 070125 in Figure [9j This is the only band affected 
by IC emission. In the upper panel of Figure [HI we plot the model assuming only synchrotron 
radiation, while the lower panel incorporates IC scattering as well. The wind model in a 
pure synchrotron fit has very unphysical parameters (Table [5} ; hence we do not consider 
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it further. In the ISM model, the observations at 4 < t < 10 d do not fit the data very 
well without any IC emission. The fit improves significantly for both environments when we 
incorporate IC effects. The IC component seems to raise the X-ray flux roughly almost the 
same time the jet break becomes visible in the optical bands. This could well explain the 
jet break at a later stage in the X-rays. We will discuss this further in §51 



5. Discussion 



With our comprehensive broadband models in hand, we now turn to some of the ques- 
tions raised in the previous sections. 



5.1. Is Inverse Compton scattering delaying the jet-break? 

Here we examine the IC scattering effect on the GRB afterglow lightcurve in the X- 
ray band. We adapt an approach in which we use only the synchrotron model for the 
GRb afterglow and derive various parameters such as E, p, e e , e#,, density etc from the 
broadband data fitting. In this approach, we force the broadband jet break to be fixed on 
the day of the optical jet break from our analytical fits i.e. on day 3.7. We then use these 
parameters to derive the light curve purely due to IC effect. For reasons noted earlier we 
confine our discussion to the ISM model. Here, we assume that the spectrum due to IC 
scattering has the same shape as that of the synchrotron model. Hence, IC spectrum in the 
X-ray band is 

\ -i/2 

pic l \ , ,IC / , , / , ,ic 

r ">^ IC / ' C ^ ^ m 1 



AC 



-P/2 /,,IC\ -1/2 ( 14 ) 



jylC I I jylC / ' m ' 



Here F^ ax is IC peak flux, z/* c is IC cooling frequency and v > v l £ is IC peak frequency. 

The best fit parameters from the synchrotron broadband model fit are: -E52 = 2.98, 9j = 
0.23 rad, p = 2.27, n = 15.7 cm" 3 , e e = 0.275, e B = 0.274, t c = 7.8d and tj = 3.7d. 
Here t c is the transition time from the fast cooling to the slow cooling state. E52 is the 
isotropic-equivalent ki netic energ at t = t c . IC scattering delays the cooling time by a 



fraction (1 + e e /e#) 2 (jSari &: Esinll200ll ). i.e. the cooling time in presence of the IC effect 
is = 3 Id. Thus, for the time span of our observations, the afterglow remains in the fast 
cooling state and we will use formulation in this regime. 
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Using the formulation described in I Sari &: Esinl (120011 ) and lWei &: Lul (119981 ). we estimate 
the time (he) when IC scattering starts becoming important at 1.486 keV (the X-ray light 
curve frequency). Using above best fit parameter values, we obtain t IC = 2.8 day. The IC 
light curve will satisfy the pre-jet break condition in the timerange of tic < t < Uet, i-e. 
between day 2.8 — 3.7. The light curve will follow post-jet break formulation from day 3.7 
onwards. 



The flux density, size, Lorentz factor and cooling frequencies derived on day 2.8 are: 
R = 2.86 x 10 17 cm, F max = 24.2 mJy, v m = 4.09 x 10 12 Hz and v c = 5.17 x 10 n Hz. 
Hence, the derived IC parameters on day 2.8 are: F^ x = 0.024 /xJy, u~ = 27^i/ m = 
2. 42 x 10 18 Hz, and v l ? = 2^lv c = 23.9 x 10 16 Hz. Here 7 m and 7 C parameters are defined 



m 



Sari &: Esinl (120011 ). The time dependences of various parameters in the pre-jet break 



epoch are: R oc t 1/4 , T oc r 3/8 , F max oc t° 



v m oc 



t~ 3 /\ v r oc rVa F™ 



or f 1 / 4 j/ IC oc 
max ^ 1 i u m ^~ 



IC 



< V < 



t 9 / 4 , and ul c oc t 1 ^ 4 : . Our frequency of observation [y = 3.59 x 10 17 Hz) satisfies v\ 

condition at t — 2.8d. Therefore, using Eq. [Mj the 1.5 keV light curve for IC scattering 



between day 2.8 — 3.7 becomes 



Fl c (t) = 0.0079 



2.8d 



1/8 



/iJy. 



(15) 



From day 3.7 onwards, we use the post-jet break formulation and derive the IC light 
curve. The time dependences of various parameters in the post-jet break regime are: R oc 
t°, T oc iT 1/2 , F max oc t~ l , u m oc t~ 2 , v c oc t°, F^ K oc t' 1 , v 1 ^ oc t~ 3 , and vf 1 oc t 1 . On 
the jet break day, we find v l £ = 1.29 x 10 18 Hz and ^ c = 3.6 x 10 16 Hz, which still satisfies 
i>l c < v < v 1 ^ condition. Therefore, we derive the light curve from day 3.7 onwards using 
Eq. Hi as 

Fl c (t) = 0.0082 l — \ /iJy. (16) 

The IC frequency u£ reaches the X-ray observation frequency of v = 3.594 x 10 17 Hz on day 
5.7. Hence, the above light curve is valid until day 5.7. 

After day 5.7, the IC flux density in the v > i/£ regime (Eq. [141) . which gives the light 
curve in this regime as follows: 

-2.4 

-.IC/ 



F^(t) = 0.0066 ( — ) /.Jy. (17) 



Combining Eqs. (fTBT) . (fTol) and (fTTl) . we plot the total IC light curve in Fig. (TUl The 
figure clearly shows that IC scattering flattens the light curve and delays the jet break to 
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later time. What happens after day 9 — 10 is unknown due to the lack of X-ray detection of 
the afterglow. 

While we have strong evidence that IC effects delay the X-ray jet break in GRB 070125, 
we would like to know if this effect could be seen in other GRBs as well. In the pre- Swift era, 
the X-ray data were not sufficiently well sampled to search for jet breaks, and so collimation 
corrections were almost exclusively calculated in the optical bands. In the Swift era, with a 
plethora of well-sampled XRT light curves, we may be missing the jet break due to IC effects 
in many GRBs. The radio afterglow is not useful in determining the jet break, since GRBs 
most likely scintillate in radio bands at such early times. This makes the optical the unique 
bandpass in which the real jet break can be determined unambiguously. 

The inverse Compton effect is most prominent in a dense medium. Our radio observa- 
tio ns have already sh own that GRB 070125 resides in a very dense medium. It has been shown 



by IWei fc Lul (120001 ) that IC is important in relativistic ejecta and even in non-relativistic 
ejecta in high density. Harrison et al. (2001) found good evidence for IC production of X-rays 
in GRB 000926 and deri ved a high circumbu rst density, 30 cm" 3 , comparable to the density 



we find for GRB 070125. ICorsi fc Pirol (120061 ) have shown that the late time flattening in the 
X-ray light curve of XRF 050406 can be explained as an effect of IC scattering. For inverse 
Compton scattering to play an important role, the electron energy density fraction (e e ) must 
be larger than the magnetic energy density fraction (es). 

Recently, ? has discussed chromatic breaks occurring due to scattering of the forward- 
shock synchrotron emission by a relativistic outflow located behind the leading blast-wave. 
This model may have X-ray jet breaks showing up at later times than the optical breaks. 
However, this model requires a long-lived central engine, which may not be the case with 
most of the GRBs. 



5.2. An Emerging Class of Hyper-Energetic (E > 10 52 erg) GRBs? 



With GRB 070125, we now have found three Swift events with t otal energy release i n 
excess of 10 52 erg: GRB 050904 jFrail et all 120061 ) and GRB050820A JCenko et al.l l2006bh . 
While both GRB 050904 and GRB 050820A appear to have exploded in an dense circum- 
burst medium, the lack of a bright radio afterglow from GRB 050820A indicates a more 
typical environment. Moreover, with the exception of the total energy release, other param- 
eters d erived from broadband modeling are in line with previous studies of less energetic 
GRBs ( IPanaitescu fc Kumar! 1200 It lYost et al.l 120031 ). It seems likely, therefore, that some 
factor intrinsic to the progenitor system is responsible for the large energy release. 
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At first blush, it seems surprising that Swift has detected three of the most energetic 
GRBs ever. With its increase high-energy sensitivity, Swift should preferentially select GRBs 
at the low end of the flu ence distribution. We not e, however, that a strong selection bias 
exists. As first noted by iKocevski fc Butler! (120071 ). in many Swift X-ray light curves, the 
last XRT measurement is not sufficient to rule out a collimation-corre cted prompt energy 
release of ~ 10 51 erg. Similarly, in the optical bandpass, iDai et al.l (120071 ) have shown that at 
least some jet breaks occur at late times beyond the sensitivity of medium aperture facilities. 

While a detailed discussion of the relative rates of hyper-energetic events is still pre- 
mature, it is clear at this point that, at the very least , the prompt 7-ray energ y distribu- 
tion is significantly broader than previously believed (IKocevski &: Butlerl 120071 ) . Coupled 
with the recent controversy surrounding the y alidity of the many high-energy correlations 



Butler et al.l 120071 : IWillingale et al.l 120071 ). we believe the future utility of GRBs as 



(e.g. 

cosmological probes is significantly lessened. 

Even more importantly, however, hyper-energetic GRBs have important co nsequences 
for pro genitor models. Sustained engine activity has been seen now in many GRBs (IBurrows et al 
2005b|). This poses a problem for the collapsar model, as the duration of the centr al engine 
shoul d not significantly exceed the accretion time scale onto the remnant black hole (jWoosley 
19931 ). Late-time engi ne activity i s naturally accommodates by models in which the central 
object is a magnetar (lUsovl Il992l ). The existence of hyper-energetic GRBs, however, is a 
direct and severe challenge to the magnetar model. 

With the current rate of hyper-energetic events (~ 1 yr _1 ), coupled with the difficulty in 
measuring late jet breaks for more typical Swift events, future prospects look grim. However, 
the impending launch of GLAST offers a new hope in the study of GRB energetics. Much 
like blazars, those GRBs capable of producing GeV photons detectable by the Large Area 
Telescope should be the most energetic and narrowly beamed events. Together, synergistic 
GLAST and Swift observations in the coming years should be able to shed light on the 
opening angles and energy release of a large sample of GRBs. 



5.3. The Unusual Environment of GRB 070125 



Our observations presented here, particularly the bright, self-absorbed radio afterglow, 
indic ate GRB 070125 ex ploded in a dense circum-burst medium. In a separate work, how- 
ever, |Q3nko_j3t_aL| (120081 ) have reported spectrosc opic observat i ons in dicating an environment 



almost completely devoid of absorbing material. ICenko et al.l (120081 ) furthermore find no ev- 
idence of an underlying host galaxy to deep (R > 25 mag) limits. We briefly reiterate here 
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the resolution of this apparent paradox. 



The key to understanding the environment of GRB 070125 is to note that the broadband 
afterglow emission and the super-posed spectroscopic absorption features derive from distinct 
physical regions. Afterglow emission is ca used by elect rons in the circum-burst medium 
accelerated by the outgoing blastwave (e.g. |Piranll2005bl ). These electrons reside relatively 
close to the explosion center, typically at radii r < 1 pc. 

We have strong evidence, however, that the absorption features seen super-posed on 
GRB afterglow spectra derive from material at significantly larger distances from the ex- 
plosion site: r ~ 1 kpc, or within the host galaxy ISM. Evidence in support of this large 
distance comes primarily from two lines of argument. First, the presence of Mg I indicates 
a large distance from the explosion site, as the first ionization energy of Mg falls below 1 
Ryd, and thus any Mg near th e GRB would be ionized to at least Mg II ( IProchaska et al. 
20071 ). Second, IVreeswijk et al.l (120071 ) have reported the detection of variability in the fine 
structure levels of Fe II from UV pumping for GRB 060418. By modeling the variability over 
time, they were able to measure the GRB-absorber distance: d = 1.7 ± 0.2 kpc. 

While this explains the apparent density paradox, we are still left to explain how a 
dense circum-burst medium could be embedded in such a tenuous ISM. Taking a clue from 
the lack of an underlying host detection, we suggested GRB 070125 may have exploded 
in a dense stellar cluster enriched by galaxy interactions. Star fo rmation in such ext reme 
environments can be seen in the local universe (e.g. Tadpole galaxy; iJarrett et al.ll2006l ). and 
under our hierarchical picture of galaxy formation, such interactions should have occurred 
more frequently at z > 1. The repor t of the detection of a faint source at the afterglow 
location at late times (IDai et al.l 120071 ) may call this interpretation into question, although 
it is unclear whether this emission is attributable to the fading afterglow or an underlying 
host (or some combination thereof). Regardless, future high-resolution imaging (i.e. HST) 
seems worthwhile to pin down the environment of this truly unique event. 



5.4. Early radio emission by the reverse shock? 

GRB 070125 was not detected at t ~ 1.5 d with the WSRT, and then was detected by 
both the VLA and the WSRT around day 5. It remained below detection level for the next 
15 days, before rebrightening on day 22. We explore the possibility that the flux from the 
GRB at t ~ 5 d could be emission from a reverse shock. We first consider the possibility 
that the non-detection was caused by modulations due to scintillation. 

The modulation in flux density due to the refractive scintillation can decrease the flux 
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density at the most up to A/„ ~ 100 fiJy. To check for this possibility, we combined all the 
eight observations taken during the 15 day non-detection phase. This vastly improved the 
signal-to-noise ratio. The flux density we obtained at the GRB position is 70 ±25 /xJy, much 
lower than scintillation can explain. 

According to the internal-external shock model for GRBs, the prompt emission is pro- 
duced by internal shocks within a relativistic outflow, while the afterglow is produced by 
external shocks with the interstellar medium. The reverse shock has a much lower temper- 
ature than that of the forward shock so it radiates at considerably lower frequencies. In 
this scenario using th e ISM model, the reverse shock emission peaks in the optical band at 
JNakar fc Piranll2005h : 

t Q = max(A/c, t dec ). 

For GRB 070125 we calculate t Q ~ 30 sec. This corresponds to the time of peak in the radio 
band to be: 

, ... <to_ 

''radio 



^radio 



Here v r a is the synchrotron self-absorption frequency, which can be written as (jNakar fc Piran 



20051 ): 

v r a (t ) = 6 x 10 12 Hz[(l + z)~( p+6 )/ 8 ef_\e^ /4 (nE 54 ) (p+6)/8 t-f +10)/8 ] 2/(p+4) - 

For values obtained from our multiwavelength analysis, we calculate v r a = 3 x 10 13 Hz and 
^radio ~ 2.25 days. Hence, the reverse shock emission peaks in the radio bands around t ~ 2 d. 
We now estimate the peak radio flux density for the reverse shock. This can be written as: 

(p-l)/2+1.3 ✓ \ (p-l)/2 
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where F Q is the peak optical flux density expressed as (jNakar fc Piranll2005l ): 

K ~ 16.6 mjy (1 + z )-^n^E^" {^j'^" (^A ^ 



For our best fit parameters, this value is ~ 1 mjy, which a gives the peak radio flux density 
to be ~ 1 /iJy, which is two orders on magnitude than the observed one. This shows that 
a reverse shock is not strong enough to explain the detection on day 5. We cannot explain 
this strange behavior. 



5.5. GRB with high radiative efficiency? 



One of the major concerns for GRB 070125 is the difference between the isotropic 7- 
ray energy obtained from the high-energy fluence and the isotropic-equivalent blastwave 
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energy obtained from our best fit model on the day when v c — u m . The isotropic-equivalent 
kinetic energy is an order of magnitude smaller than the isotropic 7-ray energy (10 54 erg, 
§3.3p . However, at very early times, the afterglow is in the fast cooling regime, where it 
undergoes significant loss of the energy because it is highly radiative . The fireball may lose 
as much as 80% of its energy during this phase ([Harrison et al.ll200ll ; I Cohen et al.lll998l ). If 
we incorporate radiative correcti ons, we may deri v e a more ac curate estimate of the actual 
isotropic blastwave energy. From IWu et al.l (120051 ); ISaril (119971 ) we find: 



for the ISM model and 



E(t) = E 



E(t) = E 



3f e 
2 



for the wind model. For the broadband modeling parameters, the isotropic kinetic energy 
one hour after the explosion is (5.04 ± 0.87) x 10 53 and (3.88 ± 3.11) x 10 55 erg for the ISM 
and the wind models, respectively. Here, the energy for the wind model is rather unphysical. 
The efficiency rj^ for the ISM model is 0.67. 



5.6. Wind model vs ISM model 

In terms of reduced \ 2 f° r best fit, wind model is slightly better. However, the wind 
model requires an electron energy density fraction close to 1, which is very unlikely. The best 
fit parameters in the wind model are rather unphysical, with less constrained boundaries. 
However, fixing the electron energy fraction to be 0.4 also gives reasonable fits. Evidence 
favoring the ISM model comes from the fireball size estimation from the scintillation data. 
However, the uncertainties in the diffraction scintillation time estimates may bring large 
uncertainties in the size estimates. Based on the energetics arguments stated above, the 
ISM model is favored over the wind model, but we can by no means definitely dismiss the 
latter. 



6. Conclusion 

GRB 070125 is one of the brightest GRBs ever detected, both in terms of its prompt 
high-energy fluence and its optical and radio afterglows. The isotropic equivalent energy for 
the GRB is 10 54 erg. This is the most extensively followed GRB in mult i wavebands in the 
Swift era. The richness of the data allowed us to derive many important properties of the 
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GRB and place useful constraints on many parameters. GRB 070125 was one of the few 
GRBs with sub-mm observations at various epochs. Our 95 GHz and 250 GHz observations 
with CARMA and IRAM, respectively, gave a robust determination of the peak flux density 
and u m , and constrained the power [yf v ) (Fig. [H]). 

Simultaneous fitting to optical and X-ray data favors a jet break at day 3.78 than the 
single powerlaw model. The evidence for the jet break is indisputable in the optical R and 
i' bands with pre-break and post-break slopes being 1.73 and 2.49 respectively. However, 
the jet break is not very prominent in the X-ray band. When we do the independent fit to 
optical and X-ray bands, the optical best fit is consistent with our joint fit. However, the jet 
break in the X-ray band is shifted to day ~ 10. Using the best fit parameters of the model, 
we show that the inverse Compton effects will dominate throughout our observations with 
pronounced effects in X-ray frequencies. These effects delay the jet break in the X-ray band. 

We had long observations of the GRB at three epochs in the 8 GHz band. These 
data gave evidence for diffractive scintillations, which gave an upper limit on the size of the 
fireball after the jet break, until the Sedov- Taylor phase started. This estimate of the fireball 
size is consistent with the one obtained from the broadband modeling in a constant density 
medium. 

We obtained synchrotron self absorption frequency estimates at various epochs from the 
VLA radio data. The evolution of the synchrotron self absorption frequency is t~°- 2A±om , 
which is consistent with the one expected (t~ ' 2 ) in the standard afterglow model. Syn- 
chrotron self absorption frequency estimates indicate that the GRB afterglow is moving in 
a dense medium. 

Our model fits could not distinguish between the ISM density profile and the wind-like 
density profile. The \ 2 fits were marginally better for the wind model but it needed an 
unphysically high electron energy fraction (~ 1). When we fixed the electron energy density 
fraction to 0.4 in the wind model, it did give decent fits and physical parameter values. 
However, the parameter values in the ISM model are more robust and change little with 
little change in the input values, unlike the wind model which is rather unstable. In both 
the ISM and the wind models, the radiative efficiency of the GRB is very high (> 60%). 

We suggest that IC scattering is a potential candidate in flattening the light curve and 
delay the jet break in other Swift events and explain the absence of a jet break in X-ray 
light curves of some of the Swift bursts. IC effects are more prominent in high density 
medium. Frequent radio measurements are necessary to measure the circumburst density of 
the medium. Hence, in absence of good radio data, one cannot determine the importance of 
IC scattering. GRB 070125 is unique because this has the richest radio data in Swift era, 
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having closely spaced X-ray light curve. 

Even though GRB 070125 has rich multiwaveband data, we could not nail down some 
of the lingering issues, such as Wind vs ISM density profile. One reason for this is that 
much of the evolution is in the jet break phase, when the ISM and Wind models have similar 
properties. Very dense samples of sensitive radio, X-ray, 7-ray and optical data from the 
very beginning until the GRB fades below the detection limit are needed. In the future, a 
combination of Swift GLAST, ALMA, EVLA, and various optical telescopes will provide 
this opportunity. 
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Table I. X-ray observations of the GRB 070125 with Swift at 3.594 x 10 17 Hz 
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10" 


12 





253 


± 





064 





73354 





081260 ± 





021280 


(3 


79 


± 





99) 


X 


10" 


12 





304 


± 





079 





73787 





071125 ± 





013240 


(3 


32 


± 





62) 


X 


10" 


12 





266 


± 





049 


1 


34534 





032554 ± 





007425 


(1 


•52 


± 





35) 


X 


10" 


12 





122 


± 





028 


1 


35747 





041983 ± 





008103 


(1 


96 


± 





38) 


X 


10" 


12 





157 


± 





030 


1 


41052 





027726 ± 





006362 


(1 


29 


± 





30) 


X 


10" 


12 





104 


± 





023 


1 


42675 





039339 ± 





009101 


(1 


84 


± 





43) 


X 


10" 


12 





147 


± 





034 


1 


47713 





032748 ± 





007624 


(1 


53 


± 





36) 


X 


10" 


12 





123 


± 





028 


1 


48861 





038279 ± 





007576 


(1 


79 


± 





35) 


X 


10" 


12 





143 


± 





028 


1 


78359 





014033 ± 





002725 


(6 


55 


± 


1 


27) 


X 


10" 


13 





053 


± 





010 


1 


89040 





013734 ± 





003671 


(6 


41 


± 


1 


71) 


X 


10" 


13 





051 


± 





014 


1 


94946 





015608 ± 





004296 


( 7 


29 


± 


2 


01) 


X 


10" 


13 





058 


± 





016 


2 


05462 





012729 ± 





002489 


(5 


94 


± 


1 


16) 


X 


10" 


13 





048 


± 





009 


2 


18650 





009233 ± 





002318 


(4 


31 


± 


1 


08) 


X 


10" 


13 





035 


± 





009 


2 


30888 





010806 ± 





002868 


(5 


04 


± 


1 


34) 


X 


10" 


13 





040 


± 





011 


2 


47372 





015085 ± 





003520 


( 7 


04 


± 


1 


64) 


X 


10" 


13 





057 


± 





013 


2 


68019 





006752 ± 





001515 


(3 


15 


± 





71) 


X 


10" 


13 





025 


± 





006 


2 


83975 





010154 ± 





002670 


(4 


74 


± 


1 


25) 


X 


io- 


13 





038 


± 





010 
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Table 1 — Continued 



Days since Counts/s (0.3-10.0 keV) Flux (0.3-10.0 keV) Flux Density 
Explosion erg s _1 cm~ 2 /zJy 



3.07971 


0.007542 ± 0.001640 


(3.52 ±0.77) x 


10" 


-13 


0.028 ±0.006 


3.38860 


0.005554 ± 0.001293 


(2.59 ±0.60) x 


10" 


-13 


0.021 ±0.005 


3.59126 


0.006728 ±0.001769 


(3.14 ±0.83) x 


10" 


-13 


0.025 ±0.007 


4.04549 


0.004496 ±0.001265 


(2.10 ±0.59) x 


10" 


-13 


0.017 ±0.006 


5.97272 


0.001798 ±0.000493 


(8.40 ±2.30) x 


10" 


-14 


0.007 ±0.002 


9.08277 


0.001052 ± 0.000293 


(4.91 ± 1.36) x 


10" 


-14 


0.004 ±0.001 


10.5438 


0.000785 ± 0.000239 


(3.67 ± 1.11) x 


10" 


-14 


0.003 ±0.001 


14.5620 


< 0.000277 


< 1.3 x 10~ 14 






< 0.001 


39.6190 


< 0.000167 


< 2.0 x 10" 15 






< 0.0002 a 



a Chandra observations, GCN 6186 
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Table 2. Optical/IR/UV observations of the GRB 070125 



Days since 


Filter 


A 


^0 


Magnitude 




Corrected 


Flux density 


Reference 9, 


Explosion 




in fim 


Jy 




mag 


mag 


fijy 




U.04uo4 


V 


n c;^ 
U.oo 


ooyu 


lo.D4 it U.UD 


fl 17 
U. 1 i 


1 Q 

lo.o 1 


1 «1 1 _l_ Q O 

101.1 it o.o 


LtV^1MOU41, OUoD 


^03^ 


B 


n A^ 

D.4:0 


A090 


1 8 Q9 4- 03 

lo.yz m u.uo 


n 9^ 

U.ZO 


lo.uy 


1 3A 3 4- 3 S. 

lo^.o in o.o 


r , PNfi03ft ftOAl 

VjV^In DUOD, DD4:1 


R7AA1 
U. U / L ± L ± l 


v 

V 


u.oo 


3^qo 
ooyu 


io 7A 4- 07 
lO. / *± HI u.u / 


n 1 7 

U. 1 ( 


1 £. ^7 
lO.O I 


1 3A 4- S ^ 
io^.u in o.o 


r , r , i N Jfi03fi ftOAi 

VjV/In DUOD, DD4:1 


fi8A33 
U.UO4:00 


Jj 


n a^ 

D.4:0 


A090 


1 Q HQ _j_ n Ofi 

ly.uo m u.uu 


n 9^ 

U.ZO 


io on 

lO.OU 


191 A 4- Q 

izi.'i m d . y 


r , Pl\Tfi03ft ftOAl 

VjrV/lNDUOD, DD4:1 


7Q1 AQ 


u 
rv 


u.uu 


3090 
ouzu 


i e _i_ n 03 
lo.oy i u.uo 


14 


1 Q AK 


lOSQlQQ 

izo.y m o.y 


l DU 


7Q37Q 
u. i yo / y 


R 


d.dd 


3090 


1 o cri _|_ n 03 

lO.dl HI U.UO 


14 


1 & 37 
lO.O I 


1 QCC p. _|_ A 3 

ioo.d m 4.0 


± DU 


7QR0Q 
u. i yuuy 


R 
rv 


u.uu 


3090 
ouzu 


1 Q ^7 _|_ n no 
lO.O i 1 u.uo 


14 


lO.^O 


1 9S 9 4- A 1 
1 Zo . z ^ L t. 1 


PRO 
1 uu 


7Q8A8 
u. i yoio 


R 
rv 


u.uu 


3090 
ouzu 


lu.Ul 1 U.UO 


14 


1 £. A7 


123 7 ± 3 7 


PfiO 
1 uu 


80087 
U.OUUO I 


R 
rv 


u.uu 


3090 

ouzu 


1 c c;o _|_ r\ 
iO.OO 1 u.uo 


14 


IO QQ 

lo.oy 


1 QQ 1 _|_ f\ p. 

ioo. i m D.D 


PfiO 
.r DU 


8099Q 
u . ouzzy 


R 


0.66 


3020 


18.60 =b — 


0.14 


18.46 


124.8 ± — 


nPNfi098 


8039R 

U.OUOZrU 


R 
rv 


u.uu 


3090 
ouzu 


lO.OO 1 u.uo 


14 


18 41 


1 30 7 4- A 9 
lOD. / m L ±.£ 


PfiO 
l DU 


8QQfi^ 
u.oyyuo 


r' 


n R3 


3fi31 


1 Q HQ _|_ f) OQ 

iy.uo i u.uy 


14 


io on 
lo.oy 


1(10 4-1 9 
iuu.y m i .z 


VjlVl wo 


n Q0833 


rv 


D.DD 


3090 

ouzu 


i c on -4- 


n i A 

U. 14 


lo.DD 


i no 7 _i_ 
1UO. / it — 


^jrOlN DU44: 


i 08QQ0 
i.uoyyu 


R 
rv 


U.DD 


3090 

ouzu 


i e 70 -1- O <\f{ 
lo. / U I U.UD 


n i a 

U. 14 


IO Kfi 
lo.OD 


I 1 Q Q _(_ 

II O.o It D.O 


PfiO 
-T DU 


1 0Q977 


R 
rv 


U.DD 


3090 
ouzu 


i q «q _i_ n 07 

lO.UO 1 U.U 1 


14 


lo.^y 


i oi 3 4- 7 A 
izi.o m / .*± 


PfiO 
-T DU 


i.uyou'i 


R 
rv 


D.DD 


3090 
ouzu 


1 Q AO _|_ f) D 1 } 

lo.uy i u.uo 


14 


IO. OO 


1 1 A S 4- 3 ^ 
1 14:. o m o.o 


PfiO 
-T DU 


1 101AR 


R 
rv 


D.DD 


ouzu 


lO.DD 1 U.U'l 


14 


IO 

lO.OZ 


118 1 4- A n 
llO.l 4. u 


PfiO 
± UU 


1.10450 


i' 


0.77 


3631 


lO.Dl 1 U.UO 


0.10 


18.51 




P60 


1 10742 




n 77 

D. 1 1 


3R31 

OUOl 


i c «i _i_ n o^i 

lO.Ul 1 u.uo 


n in 

U. 1U 


18^1 
lO.Ol 


1 A3 9 4- A 3 

14:0 . Z ^ 4. O 


PfiO 
r u U 


111 QQ7 
1 . 1 1 OO I 


., 


n 77 

D. / 1 


3fi31 

OUOl 


1 Q t;o _|_ n no 
lO.JO 1 U.UO 


10 

U. 1U 


1 £. A3 

10.4:0 


iu 94-io 

1 0*4:. Z ^ 4. O 


PfiO 
1 uu 


111 ft30 

1 . IIOOU 


1 


n 77 

D. 1 I 


3ft31 
ODOl 


lo.oy in u.uo 


10 
U. 1U 


1 ft AQ 

io.4y 


1 Afi O 4- A R 
14D.U It 4.D 


PfiO 
-T DU 


1 11929 


R 

rv 


U.DD 


3090 
ouzu 


1 s fin 4- n n°> 

lO.DU 1 U.UO 


14 


1 ^ Afi 


1 9A 8 4- 3 f\ 

1Z4.0 m o.d 


PfiO 
-T DU 


1 1 9999 


R 


U.DD 


3090 
ouzu 


io mi n O'? 

lO.DU 1 U.UO 


14 


1 S Afi 


1 9A 8 4- A 

1Z4.0 ^ 4.U 


PfiO 
l DU 


1 12516 


R 

rv 


U.DD 


3090 
ouzu 


1 S (\A 4- n Ozt 


14 


1 8 ^0 
lO.OU 


i on 9 4-/1 i 
i zu.z m 4. i 


PfiO 
-T DU 


1 12811 


rv 


U.DD 


ouzu 


lO.DU 1 U.UO 


14 


1 8 A(\ 


1 OA q 4- q a 
1Z4.0 m o.d 


PfiO 
r u U 


1 131 OS 
1 . lOlUO 


rv 


U.DO 


3090 

ouzu 


lo.DZ m u.uo 


O 1 A 
U. 14 


1 S AS 
10.40 


1 99 A 4- 3 ft 
1ZZ.4 It O.D 


PRO 
r DU 


1 13A0^ 

1 . 104UO 


i' 


n 77 

U. / I 


3ft31 
ODOl 


lo.D4 i U.U4 


10 
U. 1U 


1O.04 


1 QQ 3 4- A 8 
ioy.o It 4.o 


PRO 

-T DU 


1 1 3fiQ8 


i' 


n 77 

U. / I 


3fi31 
ODOl 


icfioin C\A 


10 
U. 1U 


1 R ^9 


1 A1 4-^1 
141 .y m o. i 


PfiO 
.r DU 


11/1 1 A O 

1.14/ 4o 


K 


U.DD 


oU2U 


1 q 'yo _l n i c 
lo. /o ± U.lo 


0. 14 


18.59 


1 1 n a _i_ i /i c 
11U.D ± 14. o 


CjOIndUoo 


1.15466 


R 


0.66 


3020 


18.63 ± 0.04 


0.14 


18.49 


121.3 ±4.6 


P60 


1.15705 


R 


0.66 


3020 


18.64 ± 0.03 


0.14 


18.50 


120.2 ±3.6 


P60 


1.16182 


R 


0.66 


3020 


18.63 ±0.03 


0.14 


18.49 


121.3 ±3.9 


P60 


1.16421 


R 


0.66 


3020 


18.66 ±0.03 


0.14 


18.52 


118.1 ±3.8 


P60 


1.16661 


i' 


0.77 


3631 


18.65 ± 0.04 


0.10 


18.55 


138.0 ±4.9 


P60 


1.16900 


i' 


0.77 


3631 


18.64 ± 0.03 


0.10 


18.54 


139.3 ±4.3 


P60 


1.17378 


i' 


0.77 


3631 


18.66 ± 0.04 


0.10 


18.56 


136.9 ± 4.8 


P60 


1.17617 


i' 


0.77 


3631 


18.68 ±0.03 


0.10 


18.58 


134.2 ±4.1 


P60 


1.17857 


R 


0.66 


3020 


18.65 ± 0.04 


0.14 


18.51 


119.2 ±3.9 


P60 


1.18096 


R 


0.66 


3020 


18.66 ± 0.04 


0.14 


18.52 


118.1 ±3.9 


P60 


1.18336 


R 


0.66 


3020 


18.64 ± 0.04 


0.14 


18.50 


120.2 ±3.9 


P60 


1.18576 


R 


0.66 


3020 


18.72 ± 0.04 


0.14 


18.58 


111.7 ±3.9 


P60 


1.18815 


R 


0.66 


3020 


18.66 ±0.05 


0.14 


18.52 


118.1 ±5.5 


P60 


1.19060 


i' 


0.77 


3631 


18.71 ± 0.04 


0.10 


18.61 


130.6 ±4.6 


P60 


1.19300 


i' 


0.77 


3631 


18.69 ±0.04 


0.10 


18.59 


133.1 ±4.5 


P60 


1.19539 


i' 


0.77 


3631 


18.67 ±0.04 


0.10 


18.57 


135.6 ±4.8 


P60 
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Table 2 — Continued 



Days since 


Filter 


A 


Fo 


Magnitude 




Corrected 


Flux density 


Reference 


Explosion 




in (im 


Jy 




mag 


mag 






1.20019 


i' 


0.77 


3631 


18.70 ± 0.04 


0.10 


18.60 


131.8 ± 4.9 


P60 


1.20261 


R 


0.66 


3020 


18.67 ± 0.05 


0.14 


18.53 


116.9 ± 5.3 


P60 


1.20500 


R. 


0.66 


3020 


18.65 ± 0.04 


0.14 


18.51 


117.2 ± 4.2 


P60 


1.20740 


R. 


0.66 


3020 


18.65 ± 0.09 


0.14 


18.51 


119.2 ± 9.7 


P60 


1.20980 


R 


0.66 


3020 


18.71 ± 0.04 


0.14 


18.57 


112.8 ± 4.0 


P60 


1.21220 


R. 


0.66 


3020 


18.72 ± 0.04 


0.14 


18.58 


111.7 ± 4.2 


P60 


1.21464 


)' 


0.77 


3631 


18.71 ± 0.05 


0.10 


18.61 


130.6 ± 5.5 


P60 


1.21704 


)' 


0.77 


3631 


18.68 ± 0.04 


0.10 


18.58 


134.2 ± 5.2 


P60 


1.21944 


1 


0.77 


3631 


18.71 ± 0.04 


0.10 


18.61 


130.6 ± 4.9 


P60 


1.22185 


1 


0.77 


3631 


18.74 ± 0.04 


0.10 


18.64 


127.0 ± 4.8 


P60 


1.22425 


1 


0.77 


3631 


18.67 ± 0.04 


0.10 


18.57 


1 O T~ 1 A O 

135.6 ± 4.8 


P60 


1.22521 


Rc 


0.66 


3020 


18.90 ± 0.2 


0.14 


18.76 


C\ A C 1 o o o 

94.6 ± 23.3 


GCN6050 


1.22910 


R 


0.66 


3020 


18.71 ± 0.05 


0.14 


18.57 


112.8 ± 5.0 


P60 


1.23151 


R. 


0.66 


3020 


18.70 ± 0.09 


0.14 


18.56 


113.8 ± 9.3 


P60 


1.23391 


R 


0.66 


3020 


18.75 ± 0.05 


0.14 


18.61 


108.7 ± 4.9 


P60 


1.23632 


R 


0.66 


3020 


18.76 ± 0.04 


0.14 


18.62 


107.7 ± 4.5 


P60 


1.27104 


Rc 


0.66 


3020 


19.00 ± 0.3 


0.14 


18.86 


86.3 ± 21.3 


GC1N6050 


1.31410 


Ic 


0.81 


2380 


18.00 ± 0.3 


0.10 


17.90 


164.7 ± 6.8 


GCN6050 


1.31410 


g' 


0.49 


3631 


19.60 ± 0.2 


0.06 


19.54 


55.4 ± 8.6 


GCN6050 


1.31410 


Rc 


0.66 


3020 


"i o o r\ 1 n 

18.80 ± 0.2 


0.14 


18.66 


1 AO O I 1 ^ O 

103.8 ± 17.8 


GCN6050 


1.35377 


V 


0.55 


3590 


19.26 ± 0.27 


0.17 


19.09 


83.0 ± 18.5 


GCN6041 


1.37590 


Rc 


0.66 


3020 


18.70 ± 0.2 


0.14 


18.56 


113.8 ± 19.5 


GCN6050 


1.37605 


R 


0.66 


3020 


19.09 ± 0.05 


0.14 


18.95 


79.6 ± 3.7 


GCN6039 


1.40437 


Rc 


0.66 


3020 


19.40 ± 0.4 


0.14 


19.26 


59.7 ± 18.8 


GCN6050 


1.79844 


R. 


0.66 


3020 


19.43 ± 0.03 


0.14 


19.29 


58.1 ± 1.6 


P60 


1.81062 


i' 


0.77 


3631 


19.49 ± 0.03 


0.10 


19.39 


63.7 ± 2.0 


P60 


1.84702 


R. 


0.66 


3020 


1 r\ ri I n no 

19.51 ± 0.03 


0.14 


19.37 


54.0 ± 1.7 


P60 


1.85907 


i' 


0.77 


3631 


19.52 ± 0.03 


0.10 


19.42 


61.9 ± 1.9 


P60 


1.90324 


R. 


0.66 


3020 


19.59 ± 0.03 


0.14 


19.45 


50.1 ± 1.3 


P60 


1.91542 


i' 


0.77 


3631 


19.58 ± 0.03 


0.10 


19.48 


58.6 ± 1.6 


P60 


1.96396 


Rc 


0.66 


3020 


19.71 ± 0.02 


0.14 


19.57 


44.9 ±0.8 


GCN6096 


1.97868 


R 


0.66 


3020 


19.62 ± 0.03 


0.14 


19.48 


48.7 ± 1.4 


P60 


1.99132 


i' 


0.77 


3631 


19.69 ± 0.03 


0.10 


19.59 


52.9 ±1.6 


P60 


2.64396 


Rc 


0.66 


3020 


20.23 ±0.1 


0.14 


20.09 


27.8 ±2.5 


GCN6047 


2.65995 


Rc 


0.66 


3020 


20.26 ± 0.11 


0.14 


20.12 


27.0 ±2.7 


GCN6047 


2.67696 


Rc 


0.66 


3020 


20.21 ± 0.11 


0.14 


20.07 


28.3 ±2.8 


GCN6047 


2.78096 


Rc 


0.64 


3020 


20.25 ± 0.11 


0.14 


20.11 


27.3 ±2.7 


GCN6047 


2.79596 


Rc 


0.66 


3020 


20.35 ± 0.12 


0.14 


20.21 


24.9 ±2.7 


GCN6047 


2.81479 


Ks 


2.22 


670 


17.86 ± 0.25 


0.02 


17.84 


49.0 ±10.1 


GCN6054 


2.81479 


J 


1.26 


1600 


18.82 ± 0.26 


0.05 


18.77 


49.7 ± 10.6 


GCN6054 


2.81479 


H 


1.66 


1024 


18.33 ± 0.25 


0.01 


18.32 


47.7 ±9.8 


GCN6054 


2.99296 


R 


0.66 


3020 


20.44 ±0.03 


0.14 


20.30 


22.9 ±0.6 


GCN6096 


3.02006 


R 


0.66 


3020 


20.44 ± 0.04 


0.14 


20.30 


22.9 ±0.8 


P60 


3.03328 


i' 


0.77 


3631 


20.47 ±0.04 


0.10 


20.37 


25.9 ±1.1 


P60 


3.62096 


Rc 


0.66 


3020 


20.80 ±0.2 


0.14 


20.66 


16.4 ±2.8 


GCN6064 
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Table 2 — Continued 



Days since 


Filter 


A 


Fo 


Magnitude 


A x 


Corrected 


Flux density 


Reference a 


Explosion 




in fim 


Jy 




mag 


mag 






3.91797 


r' 


0.63 


3631 


21.22 ± 0.09 


0.14 


21.08 


13.4 ±1.5 


CMOS 


4.03995 


R 


0.66 


3020 


21.07 ±0.07 


0.14 


20.93 


12.8 ±0.8 


GCN6096 


4.07106 


i' 


0.77 


3631 


21.03 ± 0.10 


0.10 


20.93 


15.5 ±1.0 


P60 


4.08226 


R 


0.66 


3020 


> 20.44 


0.14 


> 20.30 


< 22.5 


P60 


8.85767 


R 


0.66 


3020 


> 21.63 


0.14 


> 21.49 


< 7.5 


P60 


10.0099 


i' 


0.77 


3631 


23.75 ± 0.14 


0.10 


23.65 


1.3 ±0.2 


GMOS 


11.82110 


R 


0.66 


3020 


> 22.57 


0.14 


> 22.43 


< 3.2 


P60 


12.00096 


R 


0.66 


3020 


> 23.80 


0.14 


> 23.66 


< 1.0 


GCN6096 


21.99400 


R 


0.66 


3020 


> 25.40 


0.14 


> 25.26 


< 0.2 


LRIS 


21.99400 


g' 


0.49 


3631 


> 26.10 


0.06 


> 26.04 


< 0.2 


LRIS 


26.79396 


r 


0.67 


3631 


26.30 ±0.3 


0.14 


26.16 


0.13 ±- 


GCN6165 



Note. — P60= Palomcr 60-inch Telescope observations. 

a « GCN 6041: Swift UVOT, Marshall, F.E., vanden Berk, D.E. and Racusin, J. • GCN 6036: Swift 
UVOT, Marshall, F.E., Racusin, J. • GCN 6028: Palomcr 60in, Ccnko, S.B. and Fox, D.B. • GCN 6044: 16in 
PROMPT telescope, Haislip, J., Rcichart, D. and LaCluyzc, A. et al. • GCN 6035: TNT 0.8m telescope, 
Xing, L.P., Zhai, M., Qiu, Y.L. ct al. • GCN 6050: 50cm MITSuME Telescope, Yoshida, M., Yanagisawa, K., 
and Kawai, N., • GCN 6039: KANATA 1.5-m telescope, Uemura, M., Arai, A., and Uehara, T. • GCN 6096: 
MDM 2.4m and 1.3m telescopes, Mirabal, N., and Thorstensen, J.R. • GCN 6047: 152 cm Cassini Telescope, 
Greco, G., Terra, F., and Bartolini, C. ct al. • GCN 6054: PAIRITEL 1.3m telescope, Bloom, J.S., Starr, D., 
and Blake, C.H. • GCN 6064: 152 cm Loiano telescope, Greco, G., Terra, F., and Bartolini, C. et al. • GCN 
6165: Large Binocular Telescope, Garnavich, P., Fan, X., and Jiang, L. et al. 
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Table 3. Radio observations of the GRB 070125 



uatc oi 


Days since 


Telescope 


Frequency 


Flux density 


Error 


observation 


explosion 




in vj jtlz 


^Jy 




2007 Jan 26.82 


1.51 


WSRT 


4.86 


< 174 


87 


2007 Jan 29.32 


4.01 


VLA 


8.46 


360 


42 


2007 Jan 30.24 


4.93 


VLA 


8.46 


454 


38 


2007 Jan 30.95 


5.64 


WSRT 


4.86 


102 


26 


2007 Jan 31.07 


5.76 


VLA 


4.86 


jl41 


58 


2007 Jan 31.11 


5.80 


VLA 


8.46 


382 


52 


2007 Jan 31.76 


6.45 


GMRT 


0.61 


< 300 


150 


2007 Jan 31.83 


6.52 


MAMB02 


250 


3140 


590 


2007 Feb 01.92 


7.61 


MAMB02 


250 


1910 


720 


2007 Feb 02.06 


7.75 


VLA 


8.46 


563 


62 


2007 Feb 04.92 


10.61 


MAMB02 


250 


< 1470 


710 


2007 Feb 05.03 


10.72 


VLA 


8.46 


482 


52 


2007 Feb 05.04 


10.73 


VLA 


22.50 


1594 


70 


2007 Feb 05.07 


10.76 


VLA 


4.86 


< 132 


42 


2007 Feb 05.29 


10.98 


CARMA 


95 


2300 


700 


2007 Feb 06.24 


11.93 


VLA 


8.46 


489 


43 


2007 Feb 06.26 


11.95 


VLA 


4.86 


< 124 


34 


2007 Feb 07.11 


12.80 


VLA 


22.50 


1603 


235 


2007 Feb 07.12 


12.81 


VLA 


8.46 


405 


20 


2007 Feb 07.38 


13.07 


VLA 


14.94 


1159 


234 


2007 Feb 08.21 


13.90 


VLA 


14.94 


917 


293 


2007 Feb 08.23 


13.92 


VLA 


4.86 


< 145 


50 


2007 Feb 08.24 


13.93 


VLA 


8.46 


559 


26 


2007 Feb 08.40 


14.09 


VLA 


22.50 


1621 


151 


2007 Feb 09.21 


14.90 


VLA 


4.86 


< 108 


49 


2007 Feb 09.25 


14.94 


VLA 


8.46 


399 


69 


2007 Feb 09.27 


14.96 


VLA 


22.50 


1343 


162 


2007 Feb 09.27 


14.96 


CARMA 


95 


2300 


800 


2007 Feb 09.28 


14.97 


VLA 


14.94 


891 


129 


2007 Feb 10.37 


16.06 


VLA 


14.94 


1410 


137 


2007 Feb 10.79 


16.48 


MAMB02 


250 


2670 


930 


2007 Feb 11.08 


16.77 


VLA 


8.46 


385 


42 


2007 Feb 11.12 


16.81 


VLA 


22.50 


1367 


104 


2007 Feb 11.16 


16.85 


VLA 


4.86 


< 196 


71 


2007 Feb 12.08 


17.77 


VLA 


1.46 


< 920 


460 


2007 Feb 12.16 


17.85 


VLA 


8.46 


596 


109 


2007 Feb 12.20 


17.89 


VLA 


14.94 


1226 


155 


2007 Feb 12.21 


17.90 


CARMA 


95 


2100 


700 


2007 Feb 12.26 


17.95 


VLA 


22.50 


1222 


167 


2007 Feb 12.79 


18.48 


MAMB02 


250 


< 1270 


930 


2007 Feb 13.07 


18.76 


VLA 


1.46 


< 940 


600 


2007 Feb 13.11 


18.80 


VLA 


4.86 


< 150 


52 


2007 Feb 13.17 


18.86 


VLA 


8.46 


660 


39 


2007 Feb 13.21 


18.90 


VLA 


14.94 


1217 


197 


2007 Feb 13.25 


18.94 


VLA 


22.50 


1248 


83 
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Table 3 — Continued 



Date 01 


Days since 


Telescope 


Frequency 


Flux density 


Error 


observation 


explosion 




Ill jtlz 






2007 Feb 14.18 


19.87 


VLA 


8.46 


581 


14 


2007 Feb 16.17 


21.86 


VLA 


1.46 


< 1068 


540 


2007 Feb 16.19 


21.88 


VLA 


4.86 


< 159 


47 


2007 Feb 17.02 


22.71 


VLA 


1.46 


< 1152 


766 


2007 Feb 17.04 


22.73 


VLA 


4.86 


< 262 


64 


2007 Feb 18.16 


23.85 


VLA 


22.50 


1168 


118 


2007 Feb 18.18 


23.87 


VLA 


8.46 


303 


63 


2007 Feb 18.21 


23.90 


VLA 


14.94 


798 


182 


2007 Feb 18.27 


23.96 


CARMA 


95 


2400 


700 


2007 Feb 21.05 


26.74 


VLA 


14.94 


1101 


148 


2007 Feb 22.21 


27.90 


VLA 


4.86 


308 


78 


2007 Feb 23.07 


28.76 


VLA 


1.46 


< 984 


804 


2007 Feb 25.16 


30.85 


VLA 


22.50 


839 


73 


2007 Feb 27.20 


32.89 


VLA 


1.46 


< 978 


639 


2007 Mar 01.22 


34.91 


VLA 


4.86 


< 322 


95 


2007 Mar 01.24 


34.93 


VLA 


14.94 


432 


149 


2007 Mar 01.28 


34.97 


VLA 


8.46 


414 


66 


2007 Mar 01.30 


34.99 


VLA 


22.50 


788 


74 


2007 Mar 13.18 


46.87 


VLA 


1.46 


< 1078 


739 


2007 Mar 15.20 


48.89 


VLA 


4.86 


229 


49 


2007 Mar 15.22 


48.91 


VLA 


8.46 


443 


59 


2007 Mar 21.22 


54.91 


VLA 


4.86 


262 


52 


2007 Mar 21.26 


54.95 


VLA 


8.46 


473 


44 


2007 Mar 23.08 


56.77 


VLA 


22.50 


559 


161 


2007 Mar 25.15 


58.84 


VLA 


1.46 


< 580 


280 


2007 Apr 02.08 


66.77 


VLA 


4.86 


226 


57 


2007 Apr 02.12 


66.81 


VLA 


8.46 


345 


48 


2007 Apr 02.16 


66.85 


VLA 


1.46 


< 680 


480 


2007 Apr 02.20 


66.89 


VLA 


14.94 


530 


137 


2007 Apr 02.24 


66.93 


VLA 


22.50 


568 


137 


2007 Apr 21.13 


85.82 


VLA 


1.46 


< 1200 


570 


2007 Apr 21.18 


85.87 


VLA 


4.86 


302 


55 


2007 Apr 22.18 


86.87 


VLA 


8.46 


403 


56 


2007 Apr 22.99 


87.68 


VLA 


22.50 


615 


167 


2007 Apr 23.02 


87.71 


VLA 


14.94 


< 450 


226 


2007 May 15.03 


109.72 


VLA 


8.46 


267 


42 


2007 May 17.00 


111.69 


VLA 


1.46 


< 556 


278 


2007 May 18.81 


113.50 


VLA 


4.86 


< 290 


145 


2007 May 19.91 


114.60 


VLA 


14.94 


< 1440 


720 


2007 May 20.04 


114.73 


VLA 


22.50 


< 176 


88 


2007 Jul 04.76 


160.45 


VLA 


8.46 


145 


33 


2007 Jul 04.80 


160.49 


VLA 


4.86 


174 


45 


2007 Jul 04.85 


160.54 


VLA 


1.46 


< 162 


160 


2007 Aug 03.87 


190.56 


VLA 


22.50 


< 446 


223 


2007 Aug 11.68 


198.37 


VLA 


4.86 


< 100 


42 
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Table 3 — Continued 



Date of 


Days since 


Telescope 


Frequency 


Flux density 


Error 


observation 


explosion 




in GHz 


^Jy 


^Jy 


2007 Aug 11.72 


198.41 


VLA 


8.46 


166 


46 


2007 Aug 13.71 


200.40 


VLA 


14.94 


< 350 


169 


2007 Sep 08.64 


226.33 


VLA 


4.86 


< 141 


57 


2007 Oct 04.52 


252.21 


VLA 


8.46 


148 


60 


2007 Oct 04.60 


252.29 


VLA 


4.86 


203 


34 


2007 Nov 18.39 


297.08 


VLA 


8.46 


111 


19 


2007 Nov 18.48 


297.17 


VLA 


4.86 


161 


23 


2008 Jan 02.27 


341.96 


VLA 


8.46 


64 


18 


2008 Jan 02.35 


342.04 


VLA 


4.86 


133 


21 
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Table 4. Optical/X-ray Afterglow Model Fits 



Parameters 



Models 



ai 
a.2 
t b (days) 
X? (Total) 
X?(X - ray) 
X?(-R-band) 
X? («'— band) 



Single powerlaw Broken powcrlaw 

same slope independent slope same slope independent slope 



1.80 a 



2.23 (99 d.o.f.) 

0.88 (22 d.o.f.) 

2.22 (52 d.o.f.) 

3.43 (25 d.o.f.) 



: (1.85), a? (1.80) 



2.25 (98 d.o.f.) 

0.87 (22 d.o.f.) 

2.21 (52 d.o.f.) 

3.59 (24 d.o.f.) 



1.73 ± 0.02 
2 49+ - 85 

n Q + 2.6 

3.»_ 4 

1.30 (97 d.o.f.) 

1.32 (20 d.o.f.) 

1.54 (52 d.o.f.) 

0.77 (25 d.o.f.) 



af (1.76 ± 0.30), af (1.73 ± 0.02) 
o 40+0.86 

o o+2.6 
3.»_ 4 

1.31 (96 d.o.f.) 

1.31 (20 d.o.f.) 

1.54 (52 d.o.f.) 

0.81 (24 d.o.f.) 



Note. — All fits were performed jointly to constrain the relevant decay indices and break time to be the same 
in all bandpasses. All errors quoted are 90% confidence intervals. 

a The poor quality of the overall fit precludes meaningful estimates of the 90% confidence intervals. 
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Table 5. Comparison between Wind and ISM model 



parameters I ISM Wind Wind-fixed e e 





with IC 


no IC 


with IC 


no IC 


with IC 


no IC 


^ )iso (10 52 erg) a 


°-^°-0.24 


14.5711$ 


o.29i8£ 


u - io -0.08 


i 9 c-+0.02 
i - zo -0.17 


41 +°- 29 

u -^ 1 -0.01 


Oj (rad) 


u - zo -o.oi 


u - iZ -o.oi 


o.36±8:g 


n OO+90.0 

u -°°-0.04 


24+ - 01 

u - z ^0.02 


22+ - 01 

u - zz -0.01 


P 




o i i +0.02 

z - 1J --o.oi 


9 i 7+O.OI 
' -0.01 


2 14+0- 84 

^•-"-^-0.01 


9 9 c +0.01 
^•^-0.01 


9 10+O.OI 
z " ±o -0.03 


e e 


n 97+0.03 

u - z ' -0.01 


u - io -0.02 


0.991^ 


n 77+0.02 

U. 1 1 _ 74 


0.4 


0.4 


e B (%) 


9 77+O.44 
z -' '-0.75 


99.99l°; 


6- 7 &-7.28 


inn +o.oo 


o 70+0. 87 
°- ' °-0.03 


inn +o.oo 

luu -9.90 


n or 


42.071^ 


4 4Q+0.62 
^•^°-0.20 


2-si±8:fi 


0.9918;" 


O C9+0.02 

°-' JZ -0.43 


1 nn+0.12 

i - uu -0.02 


tjet (days) 


o fi q+o.03 




' - zy -1.09 


6-87±|:S 


q q 4 + 0.10 

^•^-o.oi 


q a 7+0.03 
°- fu -0.08 


tnonrel (days) 


pro 4 r+0.52 


yi - ' ^-4.02 


27.45liS 


30.88™ 


30.03l i 7 


36.68+ / , 2 


tcooi (days) 


8.491?;? 2 . 


4.59l°;H 


28.401°:™ 


30.89l° 8 20 


i8.i3i?:?| 


31.371};°* 


Ay 


7.4 x 10~ 8 


0.06 


0.08 


0.11 


4.0 x 10~ 7 


0.08 


Fit statistic 


592.02 


633.02 


556.26 


575.36 


584.85 


588.96 


rVd.o./. 


326.14 


411.42 


254.71 


298.52 


292.66 


335.44 


186 


186 


186 


186 


187 


186 



a Isotropic kinetic energy at the time when v c — v m . 
b Density n for ISM model in cm" 3 and A+ for wind model. 
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Fig. 1. — Plot of RHESSI data and Konus-Wind data. Their peaks are shifted by 8s because 
of the difference in light travel time between the two instruments. RHESSI recorded the 
burst explosion time (T ) to be 07:20:42 UT and Konus-Wind measured it to be 07:20:50. 
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Fig. 2. — Top: Optical (R- and i'-band) light curves of GRB 070125 as a result of joint optical 
and X-ray fits. Best-fit single power-law models are shown with dashed lines, while the broken 
power-law models are shown in solid lines. It is clear that in the optical bands, a broken power-law 
(indicating a jet break) is strongly favored. Bottom: X-ray light curve of GRB 070125, the the joint 
fit to optical and X-ray data. Again the single power-law model is shown as a dashed line, while 
the broken power-law model is shown as a solid line. Grey solid line indicates the independent fit 
to X-ray data. The independent fit is consistent with the joint fit in the optical bands but shifts 
the jet break to ~ 9 — 10 days. We discuss this in §3. II and $5j 
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Fig. 3. — Long term light curve at 8.46 GHz showing evidence of interstellar scintillation. 
This scintillation is refractive in nature and starts dominating once diffractive scintillation 
quenches (see text for more details). 



-45 - 




5000 10000 15000 
Time in seconds 



20000 




2000 4000 6000 8000 10000 12000 
Time in seconds 




5000 10000 15000 20000 25000 
Time in seconds 




2000 4000 6000 8000 10000 12000 
Time interval (AT(s)) 



0.25 




2000 4000 6000 8000 10000 
Time interval (AT(s)) 



0.25 
0.2 

0.15 
0.1 

0.05 




14 February 



5000 10000 15000 20000 
Time interval (AT(s)) 



Fig. 4. — Left panel: Lightcurves of the short timescale variations on 7, 8 and 20 February 
2007. The data were sampled at intervals of 20 mins. Right panel: Structure functions 
of the intensity fluctuations on 7, 8 and 14 February 2007. The bias of the thermal noise 
fluctuations, which adds in quadrature to the source variations in the structure function, has 
been subtracted. 
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Fig. 5. — Radio spectra of GRB 070125. We combined spectra between days 6 — 17, 18 — 34, 
56 — 86 and 160 — 200 and extracted the best fit v a . Between day 160 — 200 the spectrum 
is in the optically thin phase, v a is below the observed frequencies. The grey area indicates 
the uncertainity region for the powerlaw index slope. 
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Fig. 6. — Spectra on day 10.7 and day 23.4. We fit a smooth function for the GRB afterglow 
model and derive the -F max , v a , u m , and v c parameters. The plots at the bottom are vF v 
plots. For slow cooling the flux should peak at u m and the total power {vF v ) should peak at 
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Fig. 7. — Broadband modeling plots for the ISM model as well as the Wind model. We show 
light curves in 8.46 GHz, 4.86 GHz, R band, and i' band dataset. The darker line is for the 
ISM model and the Wind model is plotted with a thin line. The fits in the radio bands are 
not good, probably because of scintillation effects. 
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Fig. 8. — Broad band spectra at various days. The thick solid line is the ISM model and the 
thin solid line represents the Wind model. The peak flux density is well constrained due to 
submm observations. Both the models give equally good fits. 
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Fig. 9. — X-ray plot only with the pure synchrotron model. The lower figure includes inverse 
Compton effects. The IC effects start to dominate after 2.5 days. 
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Fig. 10.— Contribution of IC in the synchrotron model. The thin line represents the 
broadband model with the synchrotron component only. The thick line represents the IC 
light curve. 



